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In this article we calculate the electrical conductivity in QED using the 2PI effective action. 
We use a modified version of the usual 2PI effective action which is defined with respect to self- 
consistent solutions of the 2-point functions. We show that the green functions obtained from this 
modified effective action satisfy ward identities and that the conductivity obtained from the kubo 
, relation is gauge invariant. We work to 3-loop order in the modified 2PI effective action and show 

explicitly that the resulting expression for the conductivity contains the square of the amplitude 
f*"^ ' that corresponds to all binary collision and production processes. 
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I. INTRODUCTION 

*-£h ' Recent developments in heavy ion collisions and cosmology have stimulated interest in the theoretical understanding 
^p. of the dynamics of quantum fields out of equilibrium. The 2PI effective action is one promising method to study such 
O i' systems (for a review see [l||). 

It is well known that calculations using the 2PI effective theory involve problems with gauge invariance. It is easy to 
see how these problems arise. In general, the ward identities depend on cancellations between different topologies that 
correspond to vertex corrections and self energy corrections. In a 2PI effective theory, one uses corrected propagators 
t— ( . but not corrected vertices, and thus one does not expect the ward idenities to be satisfied. 

In this paper we study the applicability of the 2PI effective action to describe the equilibration of quantum fields. 
We consider the calculation of transport coefficients, which characterize the evolution of a system that is close to 
equilibrium over long time- and length-scales. We look specifically at the electrical conductivity of the QED plasma, 
which describes the diffusion of charge by an external electric field. Another important transport coefficient is the 
shear viscosity, which characterizes the diffusion of momentum transverse to the direction of propagation. The method 
we develop in this paper should be generalizable to the calculation of other transport coefficients. 

The study of transport coefficients has a long history. In scalar theories they have been studied using finite 
temperature quantum field theory @, H, 3, a direct ladder summation in Euclidean space Q, and 2PI effective 
action methods fg]. Gauge theories are more difficult to handle because of the subtlety of the power counting. The 
complete leading order calculation was done in Q. This calculation is not obtained directly from quantum field 
theory but is derived from kinetic theory. The equivalence of the quantum field theory and kinetic theory approaches 
has been demonstrated for scalar theories @, For gauge theories work has been done using the direct ladder 

summation in Euclidean space @, H, dynamical remormalization group methods [lfj, and 2PI methods in the large 
Nf approximation [ll|. The conductivity has recently been calculated using a diagrammatic method in which the 
ward identity is used explicitly to select contributions that will produce a gauge invariant result [121 ]. 

In this paper we demonstrate how the calculation of transport coefficients is organized in the framework of the 2PI 
effective action. We show that a gauge invariant result can be obtained by constructing a new effective theory defined 
with respect to the self-consistent solutions for the 2-point functions. This type of strategy was originally proposed 
by Baym and Kadanoff and has been discussed in the context of scalar theories in jlj]. Vertices are obtained 
by taking functional deriatives of the modified effective action with respect to the expectation values of the fields 
and the self-consistent solutions for the 2-point functions. These vertices obey a set of bethe-salpeter type equations 
which effectively restore the crossing symmetry and allow one to obtain green functions that explicitly satisfy the 
ward identity. Equivalently, when using the modified effective action to calculate the conductivity, the summation 
over ladder graphs is obtained automatically, independent of any power counting analysis. 
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In addition to discussing the general properties of the modified 2PI effective action, we perform an explicit calculation 
at 3- loop order. We obtain the integral equation that determines the conductivity. As shown in [ll|, the 2-loop term 
produces the square of the s-channel which gives the complete result at the leading-log order of accuracy. We show 
that the 3-loop term produces the missing contributions to the t- and u- channels so that the full matrix element 
corresponding to all binary scattering and production processes is obtained. We note that this is not the complete 
leading order result, since the colinear terms are not included. These terms will be present in a calculation using the 
3-loop 3PI effective action, and this work is currently in progress. 

Our calculation provides a field theoretic connection to the kinetic theory results of Q, which is useful in itself. In 
addition, it seems likely that quantum field theory provides a better framework than kinetic theory for calculations 
beyond leading order. Our results provide strong support for the use of nPI effective theories as a method to study 
the equilibration of quantum fields. 

This paper is organized as follows. 

In section Hi] we define some notation. In III A I we discuss the closed time path formalism of real time statistical field 
theory which we use throughout this paper. In III Bl we define the notation we use for propagators and vertices. 

In section [ill] we discuss the 2PI formalism. In lHI Al we give the basic structure of the 2PI effective action for QED. 
In IIIIBI we define the modified effective action as a function of the self-consistent solutions of the 2-point functions, 
and define external propagators and effective vertices. In MI CI we show that the photon propagator obeys the usual 
ward identity. In lHIDI we derive the bethe-salpeter type integral equations satisfied by the vertices. 

In section [IV] we present the calculation of the integral equation that determines the conductivity. In IIV Al we give 
the expression for the conductivity obtained from the kubo formula. In IIV Bl we obtain a bethe-salpeter type integral 
equation for the 3-point vertex from the self-consistent constraint on the 2PI effective action. In IIV CI we discuss the 
structure of each contribution to this equation and show that the complete set of diagrams includes the scattering 
amplitudes for all binary scattering and production processes. In IIV Dl we show that this integral equation is the same 
as the equation obtained in Q using kinetic theory. 

In section [V] we present our conclusions and discuss future directions. 

In Appendix [A] we give some of the technical details of the calculations presented in section IIV CI 

II. NOTATION 

A. Keldysh Representation of Real Time Finite Temperature Field Theory 

Throughout this paper we use the closed time path formulation of real time statistical field theory [IR [lfjj which 
consists of a contour with two branches: one runs from minus infinity to infinity along the real axis, the other runs back 
from infinity to minus infinity just below the real axis (for reviews see, for example, The closed time path 

contour results in a doubling of degrees of freedom. Physically, these extra contributions come from the additional 
processes that are present when the system interacts with a medium, instead of sitting in a vacuum. As a result of 
these extra degrees of freedom, n-point functions have a tensor stucture. Statistical field theory can be formulated in 
different bases, which produce different representations of these tensors. We will work in the keldysh basis. In the rest 
of this section indices in the 1-2 basis will be denoted bi and take the values 1 or 2. Keldysh indices will be written 
Cj and are assigned the values Cj = 1 := r and Cj = 2 := a. 

Throughout this section we discuss scalar fields only. The keldysh structure of the propagators and vertices for 
QED is the same as for scalars. It is straightforward to generalize the results below by including the appropriate dirac 
and lorentz structure. This will be discussed in the next section. 

We define the n-point functions in the 1-2 basis: 

G^(xi,---x n ) bl ... K :=(-i) n - 1 (V[<p{x 1 ) bx ---cf>(x n ) bn }) (1) 

The symbol V represents ordering along the closed time path. In what follows we will suppress the superscript (n) 
and the co-ordinate variables and write the n-point function as G bl ... bn , 

Vertex functions are obtained from the n-point functions by truncating external legs. In the 1-2 basis we write: 

Gbi—b n — G bl i 1 ■ ■ ■ G bnbr T 1 b " . (2) 

The vertex functions are obtained from the corresponding diagrams with an additional factor of i, This notation is 
illustrated schematically in Fig. [TJ 
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FIG. 1: Definitions of notation for propagator and vertices 



In the 1-2 basis these functions satisfy the constraints: 

E E ' ' ' E (-l) fcl+b2+ - fc " +n G blb2 ... bn = (3) 

6 1= 16 2 =1 b„=l 
2 2 2 

E E-"E r ^-6» = ° 

6 1 =16 2 =l 6 n =l 

The rotation from the 1-2 representation to the Keldysh representation is accomplished by using the transformation 
matrix: 

U K eldysh^(l-2) = ^ 1 _j ^ • (4) 

The n-point function and vertex in the Keldysh representation are given by: 

G ci - Cn =2f- 1 (7 Cl bl ---C/ c >G bl ... 6 „ (5) 

r ci-C„ = 2 f-l JJd^ . . . fJC^p&i-.&n _ 

For the 2-point green function and vertex function in momentum space we use the notation 

G ra (P) = G ret (P) ; Gar (P) = G a(fo (P) ; G rr (P) = G sym (P) = N(P) (G ret (P) - G adv (P) ) (6) 

IMP) = n ret (p) ; n ra (p) = n adv (p) ■ n aa (p) = n sym (p) = 7V(p)(n ret (p) - n adv (p)) 

where N(P) = 1 — 2n(po) and n(po) = l/(e /3p ° — 1). The expression relating the symmetric function to the retarded 
and advanced functions is the KMS condition. 

In order to simplify the notation for the vertices, we replace each combination of the indices {r, a} by a single 
numerical index. In momentum space we write: 

r ciC2 - c "( Pl , P2 , = T(i, Pl ,p 2 ,--- Pn ) (7) 

We assign the choices of the variables c\Ci ■ ■ ■ c n to the variable i using the vector 

M£)-< 2 >(::) «> 

where the symbol Cg> indicates the outer product. For each n, the zth component of the vector corresponds to a list of 
variables that is assigned the number i. To simplify the notation we drop the subscripts and write a list like r^a^ 
as rra. For clarity, the results are listed below. 

3- point functions: rrr — > 1, arr — > 2, rar — ► 3, aar — > 4, rra — > 5, ara — > 6, raa — > 7, aaa — ► 8 

4- point functions: rrrr — » 1, arrr — » 2, rarr — > 3, aarr — > 4, rrar — » 5, arar — » 6, raar — > 7, aaar — » 8, rrra — » 9, 
arra — > 10, rara — > 11, aara — » 12, rraa — » 13, araa — » 14, raaa — > 15, aaaa — > 16 



B. Propagators and Vertices 

In this section we define our notation for QED propagators and vertices. Greek letters from the beginning of the 
alphabet are dirac indices and greek letters from the end of the alphabet are lorentz indices. We use latin letters 
for 1-2 closed time path indices. We will need the fermion propagator, the photon propagator, the (two fermion - 
photon) 3-point vertex, the three photon 3-point vertex, and three different kinds of 4-point functions which couple 
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four fermions, two fermions and two photons, and four photons. Of course, the three photon vertex and all 4-point 
vertices are zero at the tree level. We will work in the high temperature limit where the masses can be set to zero. 

The concept of 'pinch singularities' plays an important role in the calculation of transport coefficients. The basic 
idea is that there is an infinite number of terms that all contribute at the same order because of the low frequency 
limit in the kubo formula f60|) . This limit produces pairs of retarded and advanced propagators which carry the same 
momenta. When integrating a term of the form J dpo G ret (P)G adv (P) , the integration contour is 'pinched' between 
poles on each side of the real axis, and the integral contains a divergence known as a 'pinch singularity' These 
divergences are regulated by using resummed propagators which account for the finite width of thermal excitations. 
This procedure introduces extra factors of the coupling in the denominators which change the power counting. As a 
consequence, there is an infinite set of graphs which contain products of pinching pairs that all need to be resummed. 

1. Propagators 

In coordinate space the propagators are written: 

S^(x 1} x 2 ); Df v {x u x 2 ) (9) 

In momentum space we will use the same expressions (to simplify the notation we do not introduce tilde's) with the 
coordinate variables {xi,x 2 ■ ■ ■ x n } replaced by the momentum variables {pi,p 2 ■ • -Pn}- Momenta are always taken 
to be incoming. Note that there is an overall delta function in momentum space so that, for example, the fermion 
propagator can be written: 

S^(P2) := S^(- P 2,p 2 ) (10) 
For the fermion propagator we use the notation: 

S~ 1 (P) = f -E(P); £(P)=7o£o(P)+pY£ s (P) (11) 

p - jM^J 

We set £ to zero expect where it is needed to regulate a pinch singularity. We write: 

S ret (P) = f Gret (P) := p2 - — (12) 

P 2 + i Sign(p )e - iSret(P) 

Sadv(P) = f G adv (P) 
G re t(P) G a dv{P) 1 



P 2 -iSign{p )e-^ adv (P) 

P(P) 



ImS(P) 



where we have defined 



p(P) := i d{P) := i {Gret{P) - G adv (P)) (13) 
:= Sign(p )A(P) 

There are no pinch terms involving photon propagators (see section IIV(1 and therefore we can write the photon 
propagator (in the feynman gauge) as: 

D r ;J(P) = -9,uG ret {P) := — (14) 

F P- + i bign(p )e 

D a ^{P) = ~9,uG adv {P) := - p2 _ / s 7 gnbo)£ 

Following the notation of [3| we write cut fermion lines and cut photon lines as follows: 

S d aP (P) = S${P) - S a a f(P) = r d(P) ; f = 6(poHP)u(P) 6(- Po )v(P)v(P) (15) 

DU P ) = D 7^ P ) D t V ( P ) = -9*» d (P) i 9»u = (-1) ftP) et\P) 
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For future use we also define the principle parts: 

Prin(P) := i (G ret (P) + G adv (P)) (16) 

S P Jt(P) = \{%$(P) + S a a f(P)) = f Prin(P) 

D p ;j n (P) = \(D r ;J(P)+D^(P)) = -^Prin(P) 

Note that to simplify the notation we use the same expressions G re t(P), G a dv{P), d(P) and Prin(P) in (fl2| . (fl4|) . 
(dU and (fl6| for fermion and photon propagators. In any equation it will be clear if a given function refers to a 
photon of fermion propagator, depending on whether the momentum variable corresponds to a fermion or a photon 
line. 

2. Vertices 

In coordinate space the vertices are written: 

A a a %(x 1 ,x 2 ,x 3 ); n a ^ T (x 1) x 2 ,x 3 ) (17) 

/ n /rSS\ab:cd / \ / n rS D\ab:cd / \ / irDD\ab;cd / \ 

\ M ) a fr lS \ x ^' x ^ x 3i x i) ; l M ) a ^ v ( x i^ x 2;x 3 ,X4); (M )^. XT {xi,x 2 ;x 3 ,x i ) 

The two terms in the first line indicate the (two fermion - photon) 3-point vertex and the three photon 3-point vertex, 
respectively. The three terms in the second line are the four fermion 4-point vertex, the (two fermion - two photon) 
4-point vertex, and the four photon 4-point vertex. The variables on each side of the semicolon indicate the legs that 
will join with a pinching pair of propagators. In section HVl we show that the pinch terms involve two kinds of vertices: 
3-point functions in which the momentum on one leg (which we take to be the middle leg) goes to zero, and 4-point 
functions which connect to two pairs of propagators with the same momenta. We write these vertices: 

KM P) ■= Hm A aw3 (j, -P - Q, Q, P) = A a ^(j, -P, 0, P) (18) 
(M ss ) ap . 7 s(j, P, K) := lim(M ss ) af3 , 7S (j, -P, P + Q, -K -Q,K)= {M ss ) a p, lS {j, -P, P, -K, K) 

with similar expressions for the vertices f2, M SD and M DD . In these expressions the index j refers to the keldysh 
component of the vertex, as defined in section III Al 

In the pinch limit, the spatial part (using [i = (0, z)) of the (2 fermion - photon) vertex can be decomposed as: 

A z (j, P) = (A(j) 7oPo + P(j)7V) P z + E(j)Y (19) 

It is easy to show that: 

% A z (j, K)% = \%A\ 3 , K) - K 2 A z (j, K) ; A z (j, K) := Tr(# A z (j, K)) (20) 

For the conductivity we will need only the quantity A z (j, K). In addition, we will show that we only need the trace 
of a pinched 4-point vertex: 

Af (P, K) := Tr(P [-iM(P, K)]l£ ) (21) 
where the quantity M(P, K) is a particular combination of keldysh components and is defined in section [IV Bl 

3. Compactified Notation 

It is important to define a simplified notation, in order to avoid a proliferation of indices which would make equations 
almost unreadable. 

In section [ill] we will use a single numerical subscript to represent all continuous and discrete indices. For example: 
a photon field is written A^(x) :— A\\ the fermion propagator is written S^g(xi, x 2 ) := 512, etc. We also use an 
einstein convention in which a repeated index implies a sum over discrete variables and an integration over space-time 
variables. 

In section |IV] we give expressions that result from summing over keldysh indices. Keldysh components and momen- 
tum variables are written explicitly. Traces are over dirac indices only. 
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III. THE 2PI FORMALISM 
A. Basic Formalism 

The partition function is defined as: 

Z[J, 77, fj, C,B}= J D[A$V] Exp i (Sci + J X A\ + fjx^x + + \c 12 A x A 2 + ^12*1*2 

Recall that using our notation, the repeated index represents a summation over all discrete indices and also a integral 
over space and time variables. In the expression above, the time integral is carried out along the closed time path 
contour so that functional derivatives produce path ordered green functions. 
The generating functional for connected fields is: 

W[J, 77, fj, C, B] = -i LnZ[ J, 77, fj, C, B] (22) 

Taking functional derivatives we obtain: 

SW A SW , SW -j 
oJx orjx or)x 

SW 1 SW 
— - = - (AxA 2 + iD 12 ) ; = i>iip2 + iSi 2 . 

0G12 2 SB 12 

where the connected propagators are defined as: 

iD 12 = (AxA 2 ) - (A 1 )(A 2 ) ; (A) = A (24) 
iSi2 = (*i*2>-<*i><*2>; (*>=V> 

The QED 2PI effective action is obtained by taking the double legendre transform of W[J, 77, fj, C, B] with respect 
to the sources. The expression is explicitly constructed so that partial derivatives with respect to the sources are zero. 
We have: 

Tty, A, S, D] = W[J, 77, fj, C, B] - J t Ax - fjx^i - - B 12 (^ 2 + iS 12 ) - ^C 12 {A X A 2 + iD 12 ) (25) 
By construction the effective action satisfies: 

(5r $r (5r 

TT = ~Jl - C\ 2 A 2 ; — — =fj\ - Bi 2 rp 2 ; = -771 + B 2 iip 2 ; (26) 

6A\ oipi Stpi 

7^12, ir^— — —i±>i 



SD 12 2 SS 12 
Eqn f25|) can be rewritten: 

T[iP,i>,A,S,D] (27) 
= S cl ^, iJ, A] + '-TthnD^ + ^Tr [{D^y 1 (D 21 - D° 21 )] - zTrLnS^ 1 - iTr [(S&- 1 (S21 - S° 21 )] + $[S, D] 

where S c i[ip, t/>, A] is the classical action and So and Dg are the free propagators given by: 

The function $>[S,D] is the sum of all 2PI diagrams. Note that for QED the function $>[S,D] is independent of 
{4>, V>, A}. The equations of motion of the mean field and the propagator are obtained from the stationarity of the 
action: 

st[^,a,s,d] _ srty,$,A, S,D] _ 

SA ' SD ~ [ ' 

5T[j>,j>,A,S,D] _ Q 6T[^,A,S,D] _ Q ST[^,j>,A,S,D] _ Q 
Sip Sip ss 
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B. The Modified Effective Action 



Practical calculations involve truncations, and it is well known that a straightforward truncation of $>[S,D] leads 
to problems with gauge invariance US] . The issue of gauge invariance can be addressed by introducing a different 
effective action defined with respect to the self-consistent solution of the propagator. We define S[ip, ip, A] and 
D[V,</U]by: 



5T[tl>,<$,A,S,D] 



s=s[i>,ii),A] 

These equations have the form of dyson equations: 



0; 



ST[ip,ip,A, S, D] 



5D 







D=D[ip,tp,A] 



$12 — (^12) 1 ~ S 12 5 S 12 : 

D12 = (D^r 1 - n 12 ; n 12 = 2 



.8$[S, D] 



SS 21 
.5$[S, D] 



S D 



(30) 



(31) 



5D 21 



S D 



Substituting the self consistent solutions we obtain the modified action: 

f [V>, 4>, A] = T[1>, A, Sty, i>, A] ,D[il>, A}} 



The equivalence of lj27j) and (|32| at the exact level was shown in [2lj . 
The external propagators are defined as 



ext \—l 



SA 2 SA 1 

We also define the following vertex functions: 



m^A]; (St?) 



8lp2 Sip; 



■f[ip,ip,A] 



A - - 

ill on 



M, 



ss 

54:21 



SA 3 ' Al32 " 
P*[S,D]. M s D 



5A 3 



.") 1:21 



o _ 1 *£>u 



(32) 



(33) 



(34) 



5-1:21 



8D128S45 



8S128S45 8S128D45 

Some useful relations can be obtained from the identities: 

S\zSz2 = 8 12 : D 13 D 3 2 = 5l2 

Differentiating (j35|) with respect to A and using f34|) gives: 

= 5ll'Ai/32'S'2'2 ; c , 12 — 2£>il'^l'32'£ ) 2'2 



M nn _ 4 8 2 *[S,D] 

m 5i;21 — 41 , 

45 



SD 12 SD 4 



SAz 



SA* 



(35) 



(36) 



Taking derivatives of the dyson equations ((3TJ) we obtain a set of bethe-salpeter type equations for the vertices. 
Using (HU and ((36j) we get: 



A 



152 



(37) 



A? 32 -< 



-SS 1 '2'SS 2 i SA 3 5D VV 8S 21 SA 3 

^■132 + * ( ^12-2'V [Sl'i&iShSw] + M\2%V [-Dl'4^435-^52' 



1 ^ 



152 



° 2 )- 1 



ni 2 ) 



2SA 3 

( 8 2 $[S,D] 5S V 2> 8 2 <S>[S,D] 6Dy 



8Sv2'SD 2 i SA 3 5D V2 ,8D 21 SA 3 J 
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These equations are represented graphically in Fig. [2l 



n> 



FIG. 2: Graphical representation of Eqn. (I37j> 



C. The External Propagator 

The external propagator can now be written in terms of these vertices. We will show below that these propaga- 
tors satisfy the usual ward identities. The basic mechanism is simple: the dyson equations IpT]) contain s-channel 
ladder resummations and the bethe-salpeter equations ([371 introduce t- and u-channels, and thus restore the crossing 
symmetry. 

We take the derivative of the modified effective action using the chain rule. In the expression below we suppress 
the arguments and write T[i/j, A, S, D] as T. We obtain: 



5 2 T 



5 2 T 



SA 2 5Ax 5D 3i SD 56 
( 8 2 T 5D : 



5D 34 SD 



<5 2 r 



s b 5A 2 SA 1 5S 34 SS 56 



SS 34 5S*56 

s b SA 2 5 Ax 



(38) 



SD 34 5Ax 



'34 



s 2 r 



SS: 



.'S4 



S 2 T 



s b 5A 2 SS 34 SD 56 



5D 34 SS. 



56 



SD 8A 2 SAi 



{1-2} 



s d 5A 2 SS 34 SAx 

Using l|27p . (|28|) . (|34| and f36|) this expression can be rewritten. We give the result separately for each term: 

S 2 T 



5 2 T 



SD 34 SD 56 
S 2 T 



SS 3 4d S56 
S 2 T SD. 



SD 34 5Ax 

s 2 r 



SSuSA-l 



s b SA 2 
$S 34 
sb8A% 



8A2SA1 
5D 34 5D 56 
s d SA 2 5 Ax 

SS 34 (5^56 

§ b SA 2 6 Ax 
{1-2} 



:S-1 



S 2 T 



SS 34 SD 5e 



6D 34 SS. 



5(i 



s b SA 2 5 Ax 



{1-2} 
{1-2} 



2i£l 4 x3.[D 33 >£l 3 ' 24 >Dv4\ + M^ 5 ^ 3 [D 55 iQ 5 ix6' D e i e ][D 33 ifl 3 i 24 i D 4 , 4 ] 
— ^413 [£33^3/24' >§4'4] — -^6^43 [$55' As/16'56'6] ["§33' A3> 24 "§4'4] 



iA° 4 x 3 [S 33 ,A 3 -24'S 4 , 4 } + {1^2} 



-^6 S 5^3^55'A 5 a6^6'6][^33^3'24^4'4] + {1^2} 



(39) 



The complicated set of terms above can be expressed in a simple form by using the bethe-salpeter equations. Substi- 
tuting lj37|) into lf39j) we arrive at the following compact result for the external propagator: 



(D^y 1 +i(A° 314 S M ,A 4 , 23 ,§ 3 , 3 ) 



From ((40]) we extract the vertex part of the 2-point function: 

nif = -»(A§ U S44'A*23'S3'3) := -iTr[A xSA 2 S] 
This result is illustrated in Fig. [31 The trace is over the indices that correspond to the closed loop. 



(40) 



(41) 



„ext 

n 12 = -i 




FIG. 3: Graphical representation of Eqn l[4T|l 



It is straightforward to show that the external propagator ([33)1 satisfies the usual ward identity: 

APS 1 )- 1 = 

We use the fact that the action and the integral measure are invariant under the transformation: 
Al — * .Ai + ftAi ; *i — ► *i - ieAi*i ; #x — ► *x + ieAx*x no sum 



(42) 



(43) 



We work to first order in A and integrate by parts where necessary so that we can extract an overall factor of A. We 
obtain: 



AZ[J, v ,f),C,B] 



(44) 



= / D[AW] 
■ Exp 



- AxftJj - ieAx(fk^i - - AiftC 12 A 2 - ieB^At - A 2 )* 1 * 2 



dx L + JxAi + fji^x + *i7/i + -C l2 AiA 2 + J Bi2*i* 2 







We rewrite this expression using the standard trick to extract the quantum fields from the path integral. We replace 
the fields by derivatives with respect to the sources, which act on the exponential factor in (|44[): 



16 IS - 18 — IS 

Ai^-—; tti-^Tjf^-; ; *i* 2 -> Tjs- 

l oJi i Orji i drji i oB\ 



(45) 



Using (|22| and lf23|) we can rewrite the derivatives with respect to the sources in terms of the expectation values of 
the fields. We obtain: 



Aift Ji - ieAi(r?iV>i - ^i»?x) - AxftCi 2 A 2 - ie-B 12 (Ax - A 2 )(V>xV>2 + i£ 12 ) = 



(46) 

We rewrite the above expression using (J26j) to replace the sources by the appropriate derivatives of the effective action. 
The result is: 

(5r ( r r \ <5r 

Ai5i— - ieAx I ^17— + 77-^1 ) + * e T7^( A i - A 2 )Si 2 = (47) 
oA x \ 8ipi dipi J 8bi2 

This result can be rewritten as a total differential by using ([43]) to obtain 

AAi = ftAx ; A-0x = -ieAx^x ; A^x = ieAx^x ; «AS' 12 = — ie(Ax - A 2 )iS' 12 ; iAD 12 = no sum (48) 

where the last two equations come directly from (f24|l . Prom l(47|) we get: 

(49) 



Ar = A Ax 77- + ADi2jjj— 



<5r (5r <5r 

AVSi-y- + A^ij— + A5i2-— = 



and thus, independent of the truncation scheme, we have: 



ST 



ST 



Using HH we obtain 



5r 

Ar = AAi-— + A^xff- + A^i— = 
8Ai oipx oipx 



„ ST . j ST ST 

^srr^rw^w^ 1 



(50) 



(51) 



We get ward identities by taking functional derivatives of lf5lj) with respect to the fields, and setting the remaining 
fields to zero. For example, differentiating with respect to A 2 and setting {A, tp, 4>} to zero gives: 



S 2 T 
SAx5A 2 



dx(Dtfy 







(52) 
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D. Bethe-Salpeter Equation 



We obtain a bethe-salpeter type integral equation for the (two fermion - photon) 3-point vertex by substituting the 
second equation in Eqn ([371 into the first and iterating. We show the result in Fig. |4] below. 



>.... . 



FIG. 4: Graphical representation of Eqn (l">3ll 



The corresponding equation can be written: 



A132 =A% sa + iMf% 4S S BB >A 6tai ,S4f4. (53) 
where we have defined the composite vertex: 

M s 12 % = Mg% 5 - l -M? 2 ? 2 , v D V6 D T ,2>M°% (54) 

To reduce the number of indices, we introduce a type of matrix notation. The indices that correspond to legs that 
are joined in a closed loop are summed over, and are not written explicitly. Using this notation Eqn (|54| is written: 

M S U % = M? 2 % - i (M SD b D M DS j (55) 

Z V / 12;45 

and Eqn (f53]) becomes: 

A132 = K% 2 + i(M ss S S (56) 

Equation (|56| is illustrated in Fig. [5l The dark box indicates the composite vertex M defined in (|55|) . For clarity, 
the indices that are not summed over are shown explicitly on the diagram. 



-- 3 + 



FIG. 5: Graphical representation of Eqn (j 56 j> 



We note that this notation does not really represent matrix multiplication because of the fact that the indices which 
are summed over cannot always be written next to each other. However, the meaning of the notation is immediately 
clear from the corresponding diagram. 

The vertex M. is obtained from (|34|) and J55J. In order to include all contributions that correspond to binary 
scattering and production processes within the 2PI formalism, we need to work to 3-loop order in the $ functional. 
The two diagrams we need are shown in Fig. [H 




FIG. 6: 2PI contributions to the $ functional 

The first graph in Fig. [6] produces the leading order contributions to the vertices and the second graph gives next-to- 
leading order terms. There are no next-to- leading order contributions to the vertices M SD and M DS . We write 



(57) 



11 



The results are shown in Fig. [7l 



-iM = 



- X + X 



FIG. 7: Contributions to the 4-point functions 

Substituting l(55l) and <[57|l into 1(56]) we iterate and keep terms that correspond to binary scattering and production 
processes. We obtain: 



A 



1.52 



= A? 32 + ( iMg s + iM^ s - Mg s S SMi Q s +\ Mg D D D M° s 



12;45 



S55' 84,14 A5'34' 



(58) 



Equation (f58|) is shown in Fig. [8l 



{c,eZ,e} 



>■■■ + x> 



+ 



c d e 



We write the final result: 



FIG. 8: Binary contributions to the bethe-salpeter equation 



A 



l:S2 



A? 3 2 + 



E 



iM. 



(i) 



See' Sa'a Ak'ju' 



12;45 J 55' °4'4 i4 -5'34 



je{a,b,c,d,e,f,g} 

the superscript j € {a, b, c, d, e, /, <?} indicates the 4-point function shown in part (J) of Fig. [U 



(59) 



IV. ELECTRICAL CONDUCTIVITY 



A. Kubo Formula 



We introduce several notational changes in this section. Keldysh components and momentum variables are written 
explicitly, and traces are over dirac indices only. All tilde's on propagators are suppressed. We also introduce a slightly 
different notation for the vertices. No factors of the coupling have been extracted from the vertices as defined in (J34j) . 
For example, using J28|) and f34|) gives A° = e7 M . In the rest of this section we extract one power of the coupling 
from both A^ and A M to obtain expressions with conventional form. To simplify the notation, we do not introduce 
additional primes on the new vertices. 

The electrical conductivity can be obtained from the kubo formula: 



qo^O 



(60) 
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p u (x,y) = {f{x)f(y)) ; f(x) = ^x)^{x) 
We can write the conductivity in terms of the polarization tensor using: 

p"(g ,0) = -4 n - t (9o,0) (61) 
The polarization tensor is given in l(4T)) . Summing over keldysh indices we obtain: 

n^(0) - *-e 2 J dPTr(A°[S s y m (P)A„(5,-P,-Q,P + Q)S ret (P + Q) (62) 
+S adv (P)A v (2, -P, -Q, P + Q)S sym (P + Q) + S adv (P) A„(6, —P, -Q, P + Q)S ret (P + Q)] ) 

Note that the surviving indices on the 3-point functions correspond to the middle legs of the 3-point vertices. Eqn. 
([62)1 can be rewritten using the kms conditions. The kms condition for the 2-point function is given in ©. Similar 
kms conditions exits for the 3-point functions. A complete list is given in The expression we need in this case is: 

r(6, Pi,p 2 ,p 3 ) = JV>(pjj)(r*(3, Pi,p 2) p 3 ) - r(2, p u p 2 , p 3 )) + iv F ( P ?)(r*(3, p u p 2 , p 3 ) - r(5, p, p 2 , p 3 )) (63) 

Substituting (J6j) and f63]) into (|62| and only including terms which contain pinch singularities gives: 

n™*(Q) = -|e 2 1 dP (iVF(Po) - N F (p + q a ))Tr(AlS adv (P)A:(3, -P, -Q, P + Q)S ret (P + Q)) (64) 

Substituting l|64p into f60]) we obtain: 

a = ~pfdP(l- n f (po))n f (p a )TT(Al ) S ret (P)ReA l (3,P)S adv (P)) (65) 

As explained in section UlAl the index '3' indicates the keldysh component that is retarded with respect to the middle 
leg. We substitute |(T2]) and ((20j) into {65]) and obtain: 

4 f . .. ReA l f3 P) 

a = -(3 / dP (1 - n / bo))n / ( Po )p(PyP 4 (P) ; P 4 (P) := OT | ' J (66) 
o j 21m2j(i J ) 

In the next three sections we obtain a self-consistent integral equation for the quantity B l (P) in (|66)l . We will show 
that this integral equation has the same form as that obtained in Q using kinetic theory. 



B. Bethe-Salpeter Equation 

From Eqn (f66|) we only need to calculate the third keldysh component of the 3-point vertex in (f59|l . We do the 
sums over keldysh indices and only include terms which contain pinch singularities. Eqn (|59| becomes: 

A QM/3 (3,P) = A° M/3 (3,P) + Yl ljdKM^(P,K)S r s lt(K)A st ^(3,K)S^(K) (67) 

j£{a,b,c,d,ej,g} 

where the 4-point function that appears in this equation is a combination of keldysh components and thermal functions: 

M (j) (P, K) = (13, P, K) + N F {K) (m^ (5, P, K) - A/ (j) (9, P, X)) 

The superscript j S {a, 6, c, (i, e, /, g} indicates the 4-point function shown in part (j) of Fig. [8l From (|66|) we need to 
find an integral equation for A l (P) = Tr(PA l (P)). We multiply both sides of (|67ll by P and take the trace. Using 
((12)1 . ((20)) and JSJ we obtain: 

A*(3,P)=Aj,(3,P) + ^ l -JdKM^{P,K)p{K)B\K) (68) 

je{a,6,c,d,e,/,g} 

Taking the real part of both sides and using the definition of B t (P) in l(66]) we obtain: 

2Im£(P)-P i (P) = ReA* ) (3,P) + ^ ^ J dK Be[M^\P,K)]p(K)^(K) (69) 

j'<E{a,fc,c,d,e,/,g} 
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C. Matrix Elements 



1. Preliminaries 

In this section we discuss the structure of each of the factors Re [ M^> (P, K) ] in (|69|) . We show that the complete 
set of diagrams produces all of the amplitudes that correspond to binary scattering and production processes. The 
calculation for each diagram is similar. We outline the procedure below. In Appendix A we give the details for two 
diagrams: M (f>) (P, K) and M (d) (P, K). 

(1) The first step is to sum over the keldysh indices. These sumations can be done by hand, but the calculation 
is extremely tedious. Instead, we use a Mathematica program. This program is described in detail in [13] and is 
available at www.brandonu.ca/physics/fugleberg/Research/Dick.html. The program can be used to calculate the 
integrand corresponding to any diagram (up to five external legs) in the keldysh, RA or 1-2 basis. The user supplies 
input in the form of lists of momenta and vertices for each propagator and vertex. 

(2) The second step is to divide the result into real and imaginary parts. The method is related to the Cutkosky 
rules at finite temperature, and is described in j23|. In Appendix A we explain in detail how the procedure works for 
two examples. 

(3) The real part of each diagram will contain two cut internal lines. Opening these two lines effectively divides each 
diagram into two pieces. The last step in the calculation is to combine the pieces from all diagrams and show that 
the total result can be written as the square of the amplitudes that correspond to binary processes. The procedure 
is as follows. For the moment, we call the momenta of the internal cut lines Pi and P2. In principle, each diagram 
M^\P,K) contains 16 terms which correspond to the 2 4 possible choices for the signs of the O-components of the 
momenta on the four on-shell lines: {p°, k , r°, r 2 }. Since P is an external variable in the integral equation (|69|) . 
we make the choice p° > 0, which leaves eight terms. Only three of these terms correspond to kinematically allowed 
binary processes, or 2 — > 2 scattering/production processes. 

To write the matrix elements in conventional form, we rewrite the four variables {P, K, P 1; P 2 } in terms of the 
new variables {P, P2, L\, L 2 } which are defined so that P + P2 = L\ + L 2 and Sign(p ) = Sign(p!]) = Sign(7j) = 
Sign(^). Using this notation, the expression for (P, K) will contain an overall factor J dL\ J dL 2 <5p+p,-Li-l 2 : = 
JdL 1 fdL 2 5 4 (P + P 2 -L 1 -L 2 ). 

For each diagram we have labeled the momenta so that the internal cut lines carry momenta {P — L, K ± L} (see 
Figs. l9lto [T6|) . In each case, there are three ways to define the variables {P, P2, Ii, L 2 } corresponding to the three 
possible ways to select P2 from the set {±K, ±(P — L), ±(K ± L)}. For the case of diagrams (b) and (d), the cut 
lines carry momenta {P — L, K — L} and the three possible choices are: 

(1) P2 = —K ; {£1, L 2 } = {P — L, —{K - L)} (70) 

(2) P 2 = K-L ; {L u L 2 } = {K, P - 1} 

(3) P 2 = -(P-I); {L u L 2 } = {K, -(K - L)} 

For diagrams for which the cut lines carry momenta {P — L, K + L} the three possible choices are: 

(1) P 2 = K- {L\, L 2 } = {P — L, K + L] (71) 

(2) P 2 = -(K + L); {Li, L 2 } = {P - L, -K} 

(3) P 2 = -(P-L); {ii, L 2 } = {K + L, —K} 

In both cases, the terms corresponding to choice (2) and choice (3) can be obtained from the term corresponding to 
choice (1) by making the changes of variables: P2 «-> —L\ and P2 — L%, We define the notation 

f dP 2 f dLi f dL 2 5 P+P2 _ Ll _ h2 ]T f(P,P 2 ;L x ,L 2 ) (72) 

J J J perms 

= j dP 2 J dL x J dL 2 5 p+ p 2 - Li -l 2 f/(P,P 2 ;Ii 1 l2) + /(P,-Ii;-P2 ) l2) + /(P,-I 2 ;Ii ) -P2)^ 



(4) Finally, we need to show that the thermal factors for each term are correct. Each term should be weighted with 
a factor that corresponds to the product of the appropriate statistical emission and absorption factors. This result is 



14 



obtained by using the identity 

1 + tf B (pS)JV>(pg) + N F (p Q 3 )N F (p° 2 ) + N B (p1)N F (p 3 ) = if p° + p° 2 + p° 3 = 
where we have used the usual definitions: 

1 



Nb(po) = 1 - 2n 6 (p ) ; n b (p ) = 
N F (po) = 1 + 2n/(p ) ; n/(po) = 



e /3Po — 1 

1 

e /3po + 1 



(73) 
(74) 



For diagrams (b) and (d) the details are given in Appendix A. The result is that the thermal factor for each diagram 
has the form of one of the two expressions below: 



Af f =n f (p )(l-n f (l 1 ))(l-n f (l° 2 )) 
A4=n / (p°)(l-n 6 (/?))(l-n 6 (/°)) 



(75) 



2. Diagram (b) 

We begin by looking at diagram (b) which is shown on the left hand side of Fig. [9l The real part of M^(P,K) 
can be written as the product of the amplitudes shown in the right hand side of Fig. [9l 



P-U 



K-L 



P-L, 



Li L .L 

X 

L 2 I., 



P-I4 



FIG. 9: The real part M (b) produces |m* + e _^ 77 



Using the first line in <[70|l we obtain (see Appendix A): 
(l-n f (p°))-Re[M b ] 

= 2 ^ I dLx JdL 2 S p+ p 2 - Li -l 2 N b {m\) ss '^ xx ' ■ (n b y s '^ xx ' A(Li) A(L 2 



(m t)-'-AA 



- (nlr'^ xx ' = e 2 <(P) (rS ret (P vl'(P 2 ) e^ix) e x J (L 2 ) 



(76) 



The amplitude (m{,) ss ~* xx corresponds to the t-channel for electron-positron production. It is easy to see that 
performing the shift of variables L\ «-> L 2 produces the w-channel. We write the result: 

(77) 



(l-n f (p°))-Re[M b ] 

^2 J dLl J dL2 S p+P2-l 1 -l 2 A4 (m^ +e _^ 



■yy ^'e+e — *77 



perms 



) A(L0 A(L 2 ) 



3. Diagram ({) 

The keldysh structure of diagram (f) is identical to that of diagram (b). In this case the fermion lines are cut 
instead of the photon lines. The diagram and the amplitudes that result from taking the real part are shown in Fig. 

Sol 



K 



P-L 



K-L 
K 



X 



L 2 



FIG. 10: The real part M (/) produces |m* +e _^ e+e - 
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Using the first line in J70|) we write the result: 

{l-n f (p°))-Re[M f ] (78) 

= 2 £ y dLi | dL 2 fr + ^_ Ll _ LlI Aff (m})"'^»« ■ (n f y s '^ s ° A(Lx) A(L 2 ) 



perms 



(m+) ss '^ S2S3 = (nj)"''^«3 = e 2 ( fi »(p)y u « a ( Ll )) D™*{P-L{) (v s *(L 2 )y»' v s ' (P 2 )) 
The amplitude (m/) ss ~* S2S3 corresponds to the t-channel for electron-positron scattering. The result is: 

(l-n / (p°))-Re[A# / ] (79) 

= 2 £ / dZi / di 2 <5 P+ p 2 _ Ll _ L2 TV/ (m^ +e -^ e+e - ■ m* +e -^ e+e -) A(ix) A(L 2 ) 

perms 

^. Diagram (d) 
Diagram (d) and its real part are shown in Fig. [TTJ 

P L K p L[ L[ p 




X fP-L. 



P+P 2 



P K ^2 L 2 L 2 P 2 

FIG. 11: The real part of M (d) gives m^ +£ _^ +e _ • m^ +e __ e+e _ 

Using the first line in ([70"1) we write the result: 

(l-n/(p°))-Re[M d ] (80) 

= ~ 2 E / di2 S p+p^-Li-l 2 Mf (ml) 88 '^ 53 ■ (n d ) ss '^ S2S3 A(L X ) A(L 2 ) 

perms 

( m t r '-^3 =e 2 ( fi .(p )7 M u «»( Ll )) D r *,(P- Ll ) (v s HL 2 )r' v s ' '(P a )) 

( nd )-'-«2«3 = e 2 (^(L 1 ) 7 I/ « S3 (i 2 )) D^(P + ft) (S S '(P 2 ) 7 !/ '^(P)) 

The amplitude (mj) ss -^ S 2 S 3 corresponds to the i-channel for electron-positron scattering and the amplitude 
(n,j) ss ^ S2S3 is the s-channel. The result is: 

(l-n/(p°))-Re[M d ] (81) 
= -2^ / di x / dL 2 5 p+ p 2 - Li -l 2 Af f (m£ e _^ e+e _ • m^ +e _^ e+e _) A(L X ) A(L 2 ) 

perms 

5. Diagram (e) 

Diagram (e) and its real part are shown in Fig. [T2l 

P P L t L[ P 

^P-K I w \ 

i -t> P-14 x V 
y\ > : _t /P+P, 



P 2 L 2 L 2 X X P 2 



FIG. 12: The real part of M (e) gives m^. e -_> e + e - • m* +e _^ e+ 
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It is easy to show that Re [M e ] = (Re [ M d ] ) * and therefore 
(l-n / (p°)).Re[M e ] 



e" 1- e — ->e T e" 



perms 



.) A(L0 A(L 2 ) 



(82) 



6. Diagram (g) 
Diagram (g) and its real part are shown in Fig. [131 



p 


L 


K 


P 


Lj 


P 


P-L 




K+L _ 




;P~Li * 


iP-Li 


P 


L 


K 


P 2 


L 2 L 2 


P 2 



FIG. 13: The real part of M (g) gives 



e e — >e e 



Using the first line in (171) we obtain 

(1 -n f (p )) - Re [M ff ] 



perms ' 



(83) 



The amplitude (m g ) ss ^ S2S3 corresponds to the i-channel for electron-electron production. It is easy to see that 
performing the shift of variables L\ <-> L 2 produces the u-channel. We write the result: 

(l-n f (p°))-Re[M 9 ] (84) 
= E I dL ^ I dL * Sp + P 2 - Ll -L 2 M f (m*t e _^ _ e _.m* e - e -^ e - e -+m^ e _^ e _ e _ -r<_ e __ e _ e _) A(L X ) A(L 2 ) 



perms ' 



7. Diagram (a) 

Now we consider diagram (a) in Fig. [68l We obtain the binary contribution to the conductivity from the diagram 
with a one loop insertion on the vertical line. This contribution and its real part are shown in Fig. [14l 



?4 



'P-K 

P-L 
^P-K 



Li L,\ 



P+0 



K+Q 



FIG. 14: The real part of M (a) gives 



Using the first line in <[70|l we write the result: 
(l-n / (p°)).Re[Ar] 

= 2 ^ jdLif dL 2 S p+ p 2 - Li -l 2 Af f (m{y 

perms 



f\Ss'— >S 2 S 3 



A(Li) A(L 2 ) 



(85) 
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The amplitude (m a ) ss ^ S2Sa corresponds to the s-channel for electron-positron scattering. We obtain: 
(1-n/O )) -Re[M a ] 

= 2 £ j dL,J dL 2 d P+ P 2 - Ll - L2 Af f (m s J +e -^ e+e --m s e+e -^ e+e -) A(L X ) A(L 2 ) 



(86) 



perms 



8. Diagram (c) 

For diagram (c) is the most difficult to handle because there are two non-zero cuts. These cuts are shown in Figs. 
[TUandQSJ We will call them the horizontal and vertical cuts. 



P-L. / K-L 



X 



L 2 



P-U 



P+Q L+0 K +Q 



FIG. 15: The horizontal cut of the real part of A/ (c) gives m*t. , ■ m" + m"'_ , in 

f o e e -t-_*-y-y e — >77 e e+— »77 e e^- >77 



I L 



P-L\ / K-L 



1 .•' \ 




P+Q 'L+Q K +Q 
1 



iP-Li 



L 2 



FIG. 16: The vertical cut of the real part of M (c) gives m'L _ _ _ • rrc" „-_>„-„- + 

1 e e — »e e e e — *e e 



. m • m - - - - 

e~ e~ — *e — e — e e — *e e 1 e — e — _,. e — e — e e — »-e e 



For the horizontal cut we use the first line in (f70jl . For the vertical cut we shift variables L — ► — L + P and use the 
first line in (f7l|) . We obtain: 

(87) 



(l-n/(p )) •Re[M c - /lorz ] 

perms ^ 

(l-n / (p°))-Re[M c -" ert ] 

£ y dLi | dL 2 6 P+P2 - Ll - L2 Af f A(£0 A(i 2 )(mLl; 



*S2S3 Si — >S 2 S 3 



+ 



perms 



(™l-HorT'^ X ' = e 2 <(P)(r S ret (P- £i) 7 1 Q , e^) e^'(L 2 ) 

(n c _, lOT . z r'^ AA ' = e 2 «(P 2 ) 7 (t^'^CP- L 2 ) 7 "') 7/3 u(P)p #(L X ) £ A '*(L 2 ) 
(mLer*)"'^' = («(P) 7 M D^iP-Lr) (u(P 2 ) 7"' u(L 2 )) 



(xic—vert ) 



ss — »AA' 



e u 



(La) 7 " «(P)) r>2*(P - L 2 ) (tt(Li) 7" u(P 2 )) 



We write the results: 

(l-n f (p )) ■Re[M c - horz ~ 
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f dL x J dL 2 8 P+ p^ Ll - L2 Mb A(Li) A(L a ) (™' f - e+ _ 



77 ^ e e+ — ^77 



perms 



(l-n f (p°))-Re[AI c ~ vert ] 

]T JdL.J dL 2 5 P+P2 _ Ll _ L2 Aff A(Li) A(L 2 ) (ro* e _^ e _ o _ • K-e-^-e- + <l e _^ 



(89) 



e~ e~ e e — >e e 



perms 



5. Combine Results 

Including the contributions from all diagrams we have 

Y Re [M^iP.K)} = Re[M e+e -_ >rr ] + Re [ M e - e -^ e - e - ] + Re [ M, 

jG{a : b,c,d,e,/,g} 

Combining ([77]), (jUJ we obtain: 

= f dLl f dL i 5 p+p-2-l 1 -l 2 Aft A(Lx) A(L 2 )| 

perms 

Combining (J84J) , ((89j) we obtain: 

(l-n / (p ))-Re[M e -e-^-e-] 

J dLl J dL ' 2 Sp+P2-l 1 -l 3 Af f A(£i) A(L 3 ) 



e+e~ — >e+ e~ 



TO e+e--^77 m e+e-^jf\ 



e e — >e e e e — >e e 



perms 



Combining ([79" ]) , fgl jl . ([H2 ]l . Ip ]) we obtain: 

(l-n / (p°))-Re[Af e+e -_ +e+e -] 

= 2 J dLl J dL2 S p+P2-l 1 -l 2 Aff A(£i) A(L 2 )|m* +e __ 



TO e+e-^e+e 



e -,- e - "" e T e -_> e -r e 



(90) 



(91) 



(92) 



(93) 



perms 



D. Integral Equation 

Finally, we substitute <(90j) , ((9TJ) , ((92j) and ((93| into the integral equation ([69]). We note that p(-P 2 )P l (-P 2 ) = 
p(P 2 )B I (P 2 ) so that both of the definitions P 2 — K and P 2 = —K produce the same overall factor. Combining all 
terms the integral equation ([69]) becomes: 

(94) 



2(1 - n f {p )) ■ ImS(P) • B l (P) 

= (l-n/(p°))-IteAj(3,P) 
1 



dP 2 / dii / d J L 2< 5 P+ p 2 _ Ll _ L2 A(P 2 )A(L 1 )A(L 2 ) £ 



perms 



\M\ B\P 2 ) 



where 



\M\ 

= Afb im:+ 



(95) 



e" 1 " e — >77 



e e — *e e 



e e — >e e 



+ 2A/> |m* +e _^ e+£ 



e^e — i-e^e 



The next step is to show that the factor M can be written in the form 



l M l 2 = 



perrns 



A/' 



/ \ ' rn e- e~ —te~ e~ m e~ e~ — >e~ er 



-Af b |m* +e _^ 77 +m^ +e -^ 77 



(96) 
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There are two points to address: 

(1) The last term in lj95j) can be written in a more symmetric way as 

2A/> K+e-^e+e- ~ "*J +e __ >e+e - | 2 -»■ TV) 

Using {H| and this result the last two terms in ([9"5| can be written: 



e^e - — ►e^e - e^e~ — >e~*-e~ 



(97) 



(98) 



perms 



The three terms contained in the sum in lj98|) correspond to the familiar results for bhabha and m0ller scattering, and 
are shown in Fig. [171 





FIG. 17: Bhabha and M0ller scattering 

(2) The first term in lj95j) is the square of the amplitude that corresponds to electron-positron production. The 
diagrams are shown in Fig. [181 



FIG. 18: electron-positron production 

The terms that give compton scattering are missing because they don't contribute to the conductivity (since they do 
not connect with the two pinched pairs of fermion propagators) . In order to write the matrix element in an symmetric 
way we must temporarily include these compton terms: 



n b (p°)(l - n,(l?))(l - n b (l° 2 ))\mt +e .^ + m« +e _ 
n 6 (p°)(l - n 6 (/?))(l - n/(l2))|m* +e _^ + m s e+e ^ 



(99) 



which allows us to write M| in the form given in Eqn. f96|) . 
Next we substitute ([96|) into ([941 and use: 



perms 



£ f(P,P 2 ,L 1 ,L 2 )- f(P,Pi,L 1 ,L 2 ) = J2 f(P,P»,LM- E f'(P,P2,L 1 ,L 2 ) 



(100) 



perms 



perms 



perms 



to obtain: 

2(1 - n/(p )) ■ ImS(P) ■ P l (P) (101) 
= (1 -n f (p ))- Re A l (3,P) + i / dP 2 J dL x J dL 2 6 P+P2 - Ll - La A(P 2 )A(L 1 )A(L 2 ) |M| 2 ^ P 4 (P 2 ) 

= (l-n / (p ))-ReA»(3,P) 

+ \J dp 2 J dL x JdL 2 |Af| 2 5 p+ p 2 _ Li _ L2 A(P 2 )A(L 1 )A(L 2 ) [B l (P 2 ] - B 1 ^) - B\L 2 
where we have used A(-X) = A(X) and P*(-X) = -B^X). 
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The last step is to rearrange (j 101 j) in the form: 

-(l-n/(p°))-ReAJ(3 J P) (102) 
= -2(1 - n/(p )) ■ Im£(P) ■ P*(P) 

+ lj dp ij dL i J dL ? |M| 2 W 2 -L 1 -L 2 A(P 2 )A(L 1 )A(i 2 ) [^(Pa)-^!)-^*^)] 

and obtain an expression for ImE(P) in terms of the squared matrix element M . Using ([3l"T) and the terms shown 
in Fig. [6] we have two contributions to S(P) which are shown in Fig. [HI 

-iZ = / > x + ■ . > ■<" > 



FIG. 19: Contributions to E(P) from 

We expand the propagators in the first diagram and extract the terms that correspond to one-loop insertions on 
each line. Combining, we obtain the three diagrams shown in Fig. [2Ql 



-i£ = 



FIG. 20: Expanded contributions to E(P) 
These three graphs can be obtained from a tadpole graph of the form shown in Fig. 




FIG. 21: Representation of E(P) in terms of a tadpole diagram 

We do not need to include all of the terms shown in Fig. |4j Using M^ b ' in Fig. [21] produces the first graph in 
Fig. [20l using produces the second graph; and produces the third graph. Summing over keldysh indices 

using the program in [121 we obtain: 

sr 4( p ) = \ E j dK ( 103 ) 

je{h,c,g} 

■(M« 75 (11,P K)S%{K) + M^. nS (7, P, K)S^(K) + M$.„ s (3, P, K)N F {K)S d Sl {K)) 

We multiply by P and take the trace. We write the fermion propagator as S(K) = If G(K) (see Eqn. (jT2j) ) and 
extract M using (plj) . Taking the imaginary part we obtain: 



hnt(P) = ~ J dKRe[M^(7,P,K)- M (j) (11, P, K) - N F (K)M (J) (3, P,K)]p(K) (104) 

Using the same method as before it is straightforward to show that the result is: 

(1 - n f (p )) • ImS(P) =-^JdP 2 I dL x J dL 2 5 P+P2 _ Ll _ L2 (105) 
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E 



777, 



,2 



e -r e — >-y~y e -r e — ^77 I 



A(P 2 )A(L 1 )A(L 2 ) 



= -\J dP * J dL i I dL ? 5 p+P2-l 1 -l 2 \M\ 2 A(P 2 )A(L 1 )A(L 2 ) 

where we have used f96|) in the last line. Substituting into (|102|) gives: 

-(l-n/(p ))-ReA* (3,P) (106) 
= lJ dP ' 2 J dLl J dL2 l M | 2 W 2 -l 1 -l 2 A(F 2 )A(L 1 )A(L 2 ) [B 4 (^) + ^(^)-S l (ii)-S i (i2)] 

This equation is exactly the same as that obtained in Q using kinetic theory. 

V. CONCLUSIONS 

In this paper we have studied the applicability of the 2PI effective action to describe the equilibration of quantum 
fields. We have calculated the QED electrical conductivity using the 2PI effective action. We have used a modified 
version of the usual 2PI effective action which is defined with respect to self-consistent solutions of the 2-point functions. 
The green functions obtained from this modified effective action satisfy ward identities and the conductivity obtained 
from the kubo relation is explicitly gauge invariant. We have shown that using this method the summation over 
ladder graphs is obtained automatically, without any power counting analysis. 

We have also done an explicit calculation at 3-loop order. We have calculated the integral equation that determines 
the conductivity and shown that the full matrix element corresponding to all binary scattering and production 
processes is obtained. Our result is complete at leading-log order but it does not contain all contributions at leading 
order since the colinear terms are not included. These terms will be present in a calculation using the 3-loop 3PI 
effective action, and this work is currently in progress. Furthermore, the method we have developed in this paper 
should be generalizable to the calculation of other transport coefficients. 

Our calculation provides a field theoretic connection to the kinetic theory results of [3], which is useful in itself. In 
addition, it seems likely that quantum field theory provides a better framework than kinetic theory for calculations 
beyond leading order. Our results provide strong support for the use of nPI effective theories as a method to study 
the equilibration of quantum fields. 



APPENDIX A 

In this appendix, we give some details of the method we use to calculate Re [ M»' (P, K) ] for each of the diagrams 
in Fig. [HI The basic strategy is explained briefly in section IIV CI Throughout this appendix, we simplify the 
notation by considering scalar fields with a cubic iteraction ((f) 3 theory). The dirac and lorentz structure of the 
fermion and photon fields will not effect the basic form of the calculation. Note that for cf) 3 theory l|2~Tj) becomes 
Re \M^'(P, K)] — ilm [Mw (P, K) ] . For simplicity we set the coupling constant to one. 

There are four basic steps to the calculation: 

Step (1): We use the mathematica program in [22] to perform the sum over keldysh indices. This program is 
described in detail in and is available at www.brandonu.ca/physics/fugleberg/Research/Dick.html. The program 
can be used to calculate the integrand corresponding to any diagram (up to five external legs) in the keldysh, RA or 
1-2 basis. The user supplies input in the form of lists of momenta, and vertices for each propagator and vertex. 

Step (2): We extract the real part of each diagram. The method is related to the Cutkosky rules at finite temperature, 
and is described in detail in [23j. The basic strategy is as follows: 

(a) We make use of the fact that it is kinematically forbidden for three on-shell lines to meet at a vertex: 

d(±(P 1 ± P 2 ))d(P 1 )d(P 2 ) = (Al) 

(b) When an integral does not contain thermal distribution functions, we can add (or subtract) terms which correspond 
to products of propagators with all poles on the same side of the real axis, since these terms are identically zero by 
contour integration. 
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We apply these rules to each diagram. The basic strategy is as follows. 
The real part is given by 

Re [ M [i) (P, K) } = [M& (P, K) + (M«) (P, K)) *} /2 (A2) 

Each diagram has four internal propagators. Terms with an even number of cut propagators are real. There are no 
terms with zero cut lines, and terms with four cut lines are zero by (|A1|) . We are left with terms with any pair of 
lines cut. Some of these terms are also zero by {33). 

To isolate the surviving terms we must take into account the fact that the momenta P and K are on-shell external 
momenta. As a consequence, Eqn (|Al|) allows us to make the replacement 

G ret (p_ R s j = Qadv _ jq = p rm (p _ jjQ ( A3 ) 

In addition we obtain relations like 

d(K - L)G ret (L) -> d(K — L)Prin(L) ; d(K — L)G adv (L) -> d(K - L)Prin(Z) (A4) 

In diagrams where all propagators carry different momenta (for example, diagram (d)), the procedure is straight- 
forward. In diagrams where more than one propagator carries a given momentum (for example, diagram (b)), one 
must be careful to show that potentially dangerous terms that contain the square of a delta function (like S(L 2 ) for 
the case of diagram (b)) do not appear. The disappearance of these unphysical terms is a well known result due to 
the KMS condition Q. 

Step (3) : We relabel the momenta so that the real part of each diagram can be written as a product of amplitudes 
that have the usual form for scattering and production amplitudes. The method is described in section IIV CI and 
illustrated below. 

Step (4) : We rewrite the thermal factors and show that they correspond to the appropriate combination of statistical 
emission and absorption factors. This is done using the identity (|73"1) . 

In the sections below, we apply this strategy to diagrams (b) and (d). 



1. Diagram (d) 

We start by looking at diagram (d) in Fig. [U First we note that it is immediately clear from looking at the diagram 
that the only possible contribution comes from the horizontal cut shown in Fig. [TTJ There is no other way to cut two 
lines without getting zero by (|A1|) . 

Step (1): We use the mathematica program in [23] to perform the sum over keldysh indices and use (| A3[) immediately. 
We separate the terms that depend on the thermal functions Nb and N F . The results are 

Re [ (P,K)} =-\ J dLd(K- L)Prin(P - K) ( A5) 

[d(P - L) (N F (fc ) - N F (fc - i )) N F ip - l ) G Tet {L) 

+ N F (fc ) N F (k - l ) G adv (P - L)G ret (L) ~ N F (fc ) N F (fc - io) G adv (L)G ret (P - L)] 

and 

Re[A~4 d) (P,K)] (A6) 
= \J dL [Prin(P - K)G adv (P - L)G let (K - L)G ret (L) + Prin(P - K)G adv (K - L)G adv (i)G ret (P - L)] 

Step (2): Using (j A4|) we can rewrite (|A5[) to obtain: 
Be[M%\p,K)] (A7) 
= ^JdLd(K- L)d(P - L)Prin( J L)Prin(P - K) (N F (fc ) N F (k - l ) - (N F (fe ) - N F (k - i„)) N F (p - l )) 

Next we rewrite (|A6|) . We use the fact that terms that contain the triple of propagators G adv (K — L)G adv (P — 
L)G Tet (L) or G Tet (K — L)G Tet (P—L)G adv (L) (and no other propagators that depend on L) can be added or subtracted 
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for free, because they give zero when the integration over Iq is done (since all poles lie on the same side of the real 
axis). Using this trick we can rewrite l|A6|) to obtain 

Re [M^ d) (P,K)] =-\J dL d(K - L)d(P - L)Prin(L) (A8) 

Combining (|A7[) and l|A8|) we obtain: 

Re[M {d) (P 7 K)] = -~ f dL d(K - L)d(P - L)Piin(L)Prin(P - K) (A9) 
• [(1 - Np (k ) (Np (k - l ) - N F (p - l )) + N F (k - l ) N F (p - l ) ] 
Step (3): Using l|13p and the first line of <f70|l . and using the notation in (|72|) . we obtain: 

Re[MW] = -\Y, f dL if dL * S p +p 2 - Li -l 2 A(Lt)A(L 2 ) (A10) 

perms 

■ Prin (P - Li) Prin (P + P 2 ) (l + N F (N F (l° 2 ) - N F (p° 2 )) - N F (l° 2 ) N F (p° 2 )) 
Step (4): Using identities of the form f73|) it is easy to show that this result can be written: 

(1 -n f (p )) -Re [M (d) ] = -2 ^ j dL x f dL 2 8p+p 2 -l 1 -l 2 Aff A (L\) A (i 2 ) Prin(i)Prin (P + P 2 ) (All) 

perms 

This equation agrees with lf81j) . except for the dirac structure. 

2. Diagram (b) 

From inspection of diagram (b) (Fig. [9} there are two possible pairs of propagators that could contribute to the 
real part: one pair corresponds to the horizontal cut and the other pair to the vertical cut. Note however that the pair 
of propagators involved in the vertical cut both carry the momentum L and thus cutting both of these propagators 
would produce a divergent term of the form S(L 2 ) 2 . We show below that the coefficient of this term is identically 
zero. 

Step (1): We sum over keldysh indices and use d(L) = G ret (L) — G adv (L) to expand all factors that depend on 
d(L). The result is: 

Re [ M (b) (P,K)]=^J dL (A12) 
(N B (k Q - l ) - N F (k )) (N B (p - l ) + N F (l )) d(K - L)d(P - L)G adv (L)G ret (L) 

- ([1 - N B (k Q - lo) (Np (k Q ) ~ N F (l Q )) - N F (k ) N F (l )} d(K - L) - G adv (K - L)) G adv (P - L)G iet (L) 2 

- ([1 - N B (k - lo) (N F (k ) - N F (l )) - N F (k ) N F (l )] d(K - L) - G ret (K - L)) G iet (P - L)G adv (L) 2 ^j 

Step (2): Using the kms condition f73|) we find that the factor in square brackets is zero. Using this result we can 
rewrite the terms proportional to G ret (L) 2 and G adv (L) 2 in the equation above as (the integrals are zero after doing 
the integration): 

\JdL (G ret (K - L)G iet (P - L)G adv (L) 2 + G adv (K - L)G adv (P ~ L)G ret (L) 2 ) = (A13) 

We rewrite the surviving terms using (from |73|) ) 

N n , N B (ko-lo)N F (ko)-l 

(l0) ~ N B (k ~lo)~Np(k ) (A14) 



24 



Re [M^(P,K)] = ~ I dL d(K - L)d(P - L) (A15) 



and obtain: 

Re [M^{P,K)] = - 

(1 + N B (p - lo) N F (k ) - N B (k - l ) (N B ( Po - l ) + N F (fc ))) G adv (L)G ret (i) 
Step (3): Using the first line in l|70p and the notation in ([721 we obtain: 

Re[M< 6 >] = -\ f dL i [ dL * Sp+p.-l ± -Lz A(ii)A(i a ) (A16) 

perms 

■ (l + N B (/§) (N B (Z°) - 7V F (pj>)) _ ^ (jO) Np ^ gadv ( p _ Li) G ret ( p _ Li) 

Step (4): Using identities of the form lj73j) it is easy to show that this result can be written: 

(1 - n/(p )) • Re [M (b) ] = -2 ^ /dLi / dL 2 N b A (Li) A (L 2 ) G adv (P - Li) G ret (P - Li) 

perms 

We compare this result with (f77| . There is an extra minus sign in (|77|) which comes from the sign associated with 
the factor v(P 2 )v(P 2 ) (see Eqn (Jl5j)). 
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